JOURNAL OF RESEARCH of the National Bureau of Standards— D. Radio Propagation 
Vol. 67D, No. 2, March- April 1963 



Height-Gain for VLF Radio Waves 

James R. Wait and Kenneth P. Spies 

Contribution from Central Radio Propagation Laboratory, 
National Bureau of Standards, Boulder, Colo. 

(Received September 24, 1962) 

The height dependence of the field strength of VLF radio waves is considered. Using 
previously developed theory, the height-gain function of the first two modes is calculated in 
terms of Airy functions of complex argument. It is indicated, for frequencies of the order of 
25 kilocycles per second, that the height-gain function reaches a maximum value at a height 
of the order of 40 kilometers when the reflecting layer is about 70 kilometers. The form of 
the height-gain function is also shown to be dependent on the finite conductivity of the 
ground. An experimental curve for 18.0 kilocycles per second based on a rocket measurement 
shows some agreement with the theory. 



1. Introduction 

Despite the extensive literature [Budden, 1962; 
Wait, 1962, etc.] on the propagation of VLF radio 
waves hardly any attention has been paid to the 
influence of varying the heights of the transmitting 
or receiving antennas. Although the theory is 
usually developed for arbitrary heights of the source 
and observer, the explicit calculations are usually 
carried out only when both heights are set equal to 
zero. Because of the long wavelength involved it is 
expected that the heights of the terminal points 
would need to be very great before any significant 
change is felt. However, with the use of high alti- 
tude rockets, field strengths can be measured right 
up to the ionosphere and, thus, the calculations are 
more than just of academic interest. 

It is the purpose of this paper to describe the results 
of some computations winch enable the height-gain 
function to be plotted for the range from the ground 
up to the reflecting layers. To simplify the calcula- 
tions the ionosphere is regarded as a sharply bounded 
medium. The excitation factors of the various 
modes are also discussed since they are closely 
related to the height-gain functions. The signifi- 
cance of the excitation factor itself has been discussed 
recently [Wait, 1962] where some explicit results 
were given and the curves of the excitation factor 
given here extend these somewhat. 

2. Formulation and Statement of Formulas 

At the risk of appearing repetitious, the formula- 
tion of the problem and the final results are briefly 
stated here. The source is a vertical electric dipole 
located on or above the surface of a smooth spherical 
earth of radius a, conductivity <r g , and dielectric 
constant e H . Spherical coordinates (/', B, 0) are 
chosen with the dipole located at r=a-\-z and 0=0. 
The ionosphere is represented as an isotropic reflect- 
ing shell located at r = a + h. The electrical proper- 



ties of this shell are not specified except to say that 
the tangential electric and magnetic fields are re- 
lated at the level r=a+h by a surface impedance Z. 
For harmonic time dependence, the radial com- 
ponent of the resultant electric field is written 
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apart from a constant factor. An expression for V 
was derived previously [Wait, 1961] as the sum of 
waveguide modes. It may be written conveniently 
in the following form 

V J^l e -"«Se""«(?»(y)(?.(y)A». ( 2 ) 
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where x=(ka/2)*d, y = (2/ka)*kzo, and y=(2/ka)* 
k{r—a). The other factors in this equation are 
discussed below. 

The complex values t n are solutions of the equation 



where 
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q i =—i(ka/2) i Zlri 0j r} Q =120w, 
q^-i(ka/2)i f-^-Y and y Q =(2/ka)*kh. 

The functions Wi(t) and w 2 (t) are Airy integral 
functions and the primes indicate a derivative with 
respect to the argument. Methods used to obtain 
numerical values of t n are described in a technical 
note [Spies and Wait, 1961]. 
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The functions G n are height-gain functions and 
they are normalized to unity for y or y equal to 
zero. Explicitly, 



G n (y)- 



where 



f(tn,V) 

"/(*., o) 



f(t n , y)=Wi(t„—y)+A(t n )w 2 (t n —y). 



(6) 
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The function A n is a modal excitation factor [Wait, 
1961]. Here it is normalized such that, in the limit 
of zero curvature and perfect ground conductivity, 
it becomes unity for all modes. Under this condition 



An=|[(f„-2 2 )- 
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In many applications at VLF the ionosphere may 
be described in terms of an effective conductivity 
parameter co r . On neglecting the terrestrial magnetic 
field and assuming that oo<^<^v where v is the collision 
frequency, it follows that u r = u>l/v where co is the 
(angular) plasma frequency. Then, to within a good 
approximation, the surface impedance is given by 
[Wait, 1962] 






('+£)©' 



being essentially independent of angle of incidence 
(or mode number) for grazing modes. 

3. Excitation Factor 

The excitation factor A n determines the efficiency 
of launching of a given mode in the earth ionosphere 
waveguide. Using (8) and the known values of 
t n [Spies and Wait, 1961], calculations of A n were 
carried out. The resulting curves of A x and A 2 are 
shown in figures la and lb for o>=<» corresponding 
to perfectly conducting ground. To simplify matters 
even further, q. t is taken to be oo which corresponds 
to a perfectly reflecting ionosphere with a local 
reflection coefficient of —1. It is immediately 
evident from these curves that the excitation of the 
least attenuated mode becomes very weak at the 
higher frequencies. However, for the range of 
frequencies and heights shown, the second mode is 
being well excited. In fact, there may be a 2 db 
increase over the excitation amplitude for a flat 
earth. 

For the cases shown in figures la and lb, A n is 
real. Actually, finite ground conductivity or losses 
in the ionospheric reflecting layer will modify A n 
slightly and cause it to become complex. To il- 
lustrate this phenomenon, Ai is calculated for various 
values of <r g for an w r value of 2X10 5 which is 
considered to be typical for daytime conditions. 
The results are shown in figures 2 a and 2b. It is 
immediately evident, on comparing the curve for 
(r g =°° with the corresponding curve for h=70 km 
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Figuees la and lb. The excitation factor for the first and 
second order modes for an earth-ionosphere waveguide with 
perfectly reflecting boundaries. 

In this case A n is real. 

in figure la, that the influence of the finite value of 
o) r is quite small. However, the effect of finite 
ground conductivity is quite noticeable. It is 
rather interesting; to see that the excitation factor 
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Figures 2a and 2b. The (complex) excitation factor for the 
first-order mode for imperfectly reflecting bou7idaries. 

is actually increased as the ground conductivity is 
decreased. This fact is not incompatible with the 
behavior of the height-gain functions as discussed 
below. 



4. Height-Gain Functions 

The height-gain factors G n (y) and G n (y) have the 
same functional form and, therefore, only one of these 
quantities need be considered. The factors G n {y) 
may be computed either directly from (6) or use may 
be made of power series developments [Wait, 1962]. 
For the results given here, the former method was 
employed. 

The results are shown in figures 3a to 6b with all 
relevant parameters being indicated on the curves. 
In figures 3a, 3b, 4a, and 4b, G n (y) is real since both 
the ground and the ionosphere are perfectly reflect- 
ing. However, when either co r or a g is finite, G n (y) is 
complex as indicated in figures 5a ? 5b, 6a, and 6b. 

The sets of curves in figures 3a and 31) are for a 
fixed height h of 70 km and, because £2— °° , the value 
of G n (y) must go to zero at 2=70 km for all the fre- 
quencies. A similar behavior is evident for figures 
4a and 4b but here the ionosphere heights are taken 
to be 60, 70, SO, 90, and 100 km in turn. When oo r 
is finite it can be seen from figures 5 and 6 that G n (y) 
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Figures 3a and 3b. The height-gain functions for the -first and 
second order modes for an earth-ionosphere waveguide with 
perfectly reflecting boundaries. 

In this case G n (y) is real. In these curves the reflecting height is fixed at 70 km. 
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Figures 4a and 4b. The height-gain functions for the first and 
second order modes for an earth-ionosphere waveguide with 
perfectly reflecting boundary. 

In these curves the frequency is fixed at 20 kc/s. 



is finite at the upper boundary and the phase is of the 
order of —90°. Such a behavior is compatible with 
the boundary condition at the sharply bounded ion- 
ized layer. 

It is seen'from the curves (e.g., figures 3a, 4a, 5a, and 
6a) that the field associated with the first order mode 
rises to a maximum at a level below the reflecting 
layer. In this sense the mode is analogous to Lord 
Rayleigh's whispering gallery (acoustical) mode 
which he observed on the inside of the dome of St. 
Paul's cathedral [Rayleigh, 1910; Budden, 1962]. 
When dealing with the earth-ionosphere waveguide 
it is appropriate to call this the earth-detached mode 
since it depends only very slightly on ground con- 
ductivity [Wait, 1962; Spies and Wait, 1961]. It 
may be observed that the form of the height variation 
is generally quite insensitive to the electrical charac- 
teristics of the reflecting boundaries. 

The effect of finite ground conductivity is seen by 
comparing figures 5a and 5b for 0>=°° with figures 
6a and 6b for <r g =5X10~ 3 . It is evident that the 
finite conductivity of the lower boundary modifies 
the amplitude of G n (y) only slightly although the 
phase is noticeably affected. This behavior is con- 
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Figures 5a and 5b. The amplitude and phase of the height- 
gain function for an imperfectly reflecting upper layer and a 
perfectly conducting lower boundary. 

The dashed curve in (a) is based on a rocket measurement at 18.0 kc/s during 
daytime for a path over sea water. 



sistent with the series development for height-gain 
functions given previously [Wait, 1962]. In par- 
ticular, it was indicated that, for low heights, the 
finite conductivity tends to reduce \G n (y)\ below 
unity and the phase of G n {y) is positive by a small 
amount. For propagation over sea water it is 
permissible to regard <r g as being infinite. 

As mentioned in the introduction, the height-gain 
functions should describe the variation of the field 
measured in a rocket as it proceeds from the ground 
up to the ionosphere. Such an experiment has 
actually been carried out by Lomax 11961] at Eglin 
Air Force Base in Florida. He employed the 18.0 
kc/s transmission from station NBA in the Canal 
Zone, Panama. At this distance, in the daytime, 
it is expected the first-order mode predominates. 
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Figures 6a and 6b. The amplitude and phase of the height- 
gain function for an imperfectly reflecting upper layer and a 
finitely conducting lower layer (typical of propagation over 
land) . 



A sketch of Lomax's results for the vertical electrical 
field is shown in figure 5a. It follows the theoretical 
curve for 15 kc/s quite closely. Over the range of 
heights from 2.5 to 60 km Lomax estimates that the 
signal strength as plotted is accurate to within 2 db. 
The slight increase at very low heights has been 
attributed to the exhaust plume which has the effect 
of increasing the effective length of the receiving 
antenna. At great heights there were also some 
departures between the calculated and experimental 
curves, but here the experimental data is not con- 
sidered to be reliable. 



5. Concluding Remarks 

The curves given in this paper, although based on 
an idealized model, should provide preliminary infor- 
mation concerning the behavior of the VLF field at 
large heights in the atmosphere. Contrary to com- 
monly accepted ideas the field may increase by a 
significant amount as the transmitting or receiving- 
antennas are raised. The effect should be par- 
ticularly noticeable at the upper end of the VLF 
band at long ranges. 

It might be mentioned that the results in this 
paper, for the height-gain functions, refer specifically 
to the vertical electric field. However, to a very 
good approximation they also hold for the horizontal 
magnetic field component. On the other hand, the 
horizontal electric field component which was not 
discussed in this paper actually is proportional to 
the derivative with respect to height of the function 
G n {y). Thus, we can expect the horizontal electric 
field to be very small on the ground and at low 
heights. However, at heights above 30 km the 
horizontal field may be quite large. 
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